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Abstract 

PHENIX and STAR data on the intercept parameter of the two-pion Bose-Einstein correlation 
functions in -^/snn = 200 GeV Au+Au collisions are analyzed in terms of different models of 
hadronic multiplicities. To describe this combined PHENIX and STAR dataset, an in-medium rf 
mass reduction of at least 200 MeV is needed, at the 99.9 % confidence level in the considered 
model class. Such a significant rf mass modification may indicate the restoration of the ?7a(1) 
symmetry in a hot and dense hadronic matter and the return of the 9th "prodigal" Goldstone 
boson. 
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Although the quark model exhibits a £7(3) chiral symmetry in the limit of massless 
up, down and strange quarks, and in principle 9 massless Goldstone modes are expected 
to appear when this symmetry is broken, only 8 light pseudoscalar mesons are observed 
experimentally. This puzzling mystery is resolved by the Adler-Bell-Jackiw Ua(1) anomaly: 
instantons tunneling between topologically different QCD vacuum states explicitely break 
the Ua(X) part of the U(3) symmetry. Thus the 9th Goldstone boson is expected to be 
massive, and is associated with the rj' meson, which has a mass of 958 MeV, approximately 
twice that of the other pseudoscalar mesons. In high energy heavy ion collisions, where a hot 
and dense medium is created, the £7a(1) symmetry of the strong interactions may temporarily 
be restored jl-3], even below the critical temperature for quark deconfmement jj]. Thus the 
mass of the r/(958) mesons may be reduced to its quark model value of about 500 MeV, 
corresponding to the return of the "prodigal" 9th Goldstone boson Here we report on an 
indirect observation of such an in-medium rf mass modification based on a detailed analysis 
of PHENIX and STAR charged pion Bose-Einstein correlation (BEC) data jsl, 6]. 

The abundance of the rj mesons with reduced mass may be increased at low p-r , by more 
than a factor of 10. One should emphasize that the rj (and rj) mesons almost always decay af- 
ter the surrounding hadronic matter has frozen out, due to their small annihilation and scat- 
tering cross sections, and their decay times that are much longer than the characteristic 5-10 
fm/c decoupling times of the fireball created in high energy heavy ion collisions. Therefore 
one cannot expect a direct observation of the mass shift of the rf (or if) mesons: all detection 
possibilities of their in-medium mass modification have to rely on their enhanced production. 
An enhancement of low transverse momentum rj' mesons contributes to an enhanced produc- 
tion of soft charged pions mainly through the v( — > i] + 7r + + 7r _ — > (tt + + 7r° + 7r _ ) + 7r + + 7r~ 
decay chain and also through other, less prominent channels. As the r/' decays far away 
from the fireball, the enhanced production of pions in the corresponding halo region will 
reduce the strength of the Bose-Einstein correlation between soft charged pions. The trans- 



verse mass (mx = ^/m 2 + p T ) dependence of the extrapolated intercept parameter A* of the 
charged pion Bose-Einstein correlations was shown to be an observable that is sensitive to 



A*(mx) data at low transverse 



such an enhanced 7/ multiplicity [7] . The proposed decrease o E 
mass has been observed both by PHENIX [5] and STAR 0, §J]. 

Our main analysis tool was a Monte-Carlo simulation of the transverse mass dependence 
of the long lived resonance multiplicities including the possibility of an enhanced rf pro- 



duction at low transverse momentum, due to a partial in-medium C/a(1) restoration and a 
related 7/ mass modification. This model and the related reduction of the effective intercept 
parameter of the two-pion Bose-Einstein correlation function was proposed first in ref. 7J 
and detailed recently in refs. 

In thermal models, the production cross sections of the light mesons are exponentially 
suppressed by the mass. Hence one expects about two orders of magnitude less rj' mesons 
from the freeze-out than pions. This suppression, however, may be moderated as a con- 
sequence of a possible rj mass reduction, and the rj' mesons may show up in an enhanced 
number. The number of in-medium rf mesons is calculated with an improved Hagedorn 
formula yielding the following 7/ enhancement factor: 

frf = — e Tc ° nd " 1 

This formula includes a prefactor with an expansion dynamics dependent exponent a ~ 
l — d/2 for an expansion in d effective dimensions 11]. As a default value, a = was taken Q 



and, for the systematic investigations, this parameter was varied between —0.5 < a < 0.5. 
Other model parameters and their investigated ranges are described as follows: T con( j in the 
above formula corresponds to the temperature of the medium when the in-medium modified 



7/ mesons are formed; its default value was taken to be T con d = 177 MeV |7] and varied 
systematically between 140 and 220 MeV. Resonances with different masses were simulated 
with a mass dependent slope parameter T e s = Tpo + fn(uT) 2 , where the default values of 



Tpo = 177 MeV and (ut) = 0.48 [12| were utilized and systematically varied in the range of 
100 MeV < Tpo < 177 MeV and 0.40 < (u T ) < 0.60 . 

Once produced, the rj' is expected to be decoupled from other hadronic matter, since its 
annihilation and scattering cross sections are very small |2j. If the 7/ mass is reduced in 
the medium, the observed 77' spectrum will consist of two components. If the of the 77' 
is large enough, it can get on-shell and escape. This will produce a thermal component of 
the spectrum. Energy conservation at mid-rapidity implies m*, 2 +Px ri ' 2 — rfir] 2 +Pt,t]' 2 - 
(In the latter equation the quantities marked with an asterisk denote the properties of the 
in-medium 77', while the ones without an asterisk refer to the free 77'.) On the other hand, 



r/'-s with p^ v , < y m rf 2 — m v ' 2 will not be able to leave the hot and dense region through 
thermal fluctuation since they cannot compensate for the missing mass {2, 3], and thus will 
be trapped in the hot and dense region until it disappears. As the energy density of the 



medium is dissolved, the effect of QCD instantons increases and the trapped rf mesons 
regain their free mass and appear at low px- 

Previous simulations [7] considered the trapped rf mesons to leave the dissolving medium 
with a negligible pr, which resulted in a steep hole in the extrapolated intercept parameter 
A*(mx) at a characteristic transverse mass of < 250 MeV {?], Q|. In this simplified 
scenario the only free parameter was the in-medium rf mass, determining the depth of the 
observed hole. In the present analysis the r/'-s from the decaying condensate are given a 
random transverse momentum, following Maxwell-Boltzmann statistics with a newly intro- 
duced inverse slope parameter -B -1 , which is necessary to obtain a quality description of the 
width and the slope of the A*(?tit) data of PHENIX and STAR in the mx ~ 300 MeV region. 
Physically, B~ l is limited by Tpo, so the trapped r/'-s may gain only moderate transverse 
momenta. Hence the enhancement will mostly appear at low px [1—3] . 

Considering the A*(mx)/A™ ax distribution, shown in Fig. [TJ a dip is observable in mea- 
sured RHIC -y/SNN = 200 GeV central Au+Au collision data at low mx values. (The A* 



values used in this analysis and their total errors are discussed in details in ref. [10j. Here 
^max j g A t (wiT) value taken at = 0.7 GeV, with the exception of the STAR data, 
where the data point at the highest mx = 0.55 GeV is considered. Note that the mx 
dependency of the A t (mx) measurements in the 0.5-0.7 GeV region is very weak.) 

We have investigated a broad class of models of resonance production, including three 
different models that produce resonances without assuming local thermalization: AMPT, 
FRITIOF and UrQMD. AMPT, known to be fairly successful in describing the HBT radii 



without assuming any in- medium mass modification [14j . is not able to describe the dataset 
shown in Fie. [Tl The trend seen in AMPT can probably be attributed to a lower effective 

n 

(ut) of the high mass halo resonances [7[. (Note that A, the Gaussian approximation of the 
intercept parameter reported by AMPT, is distinguished here from the extrapolated inter- 
cept A*, which also includes the uncertainty from the fitting form 10].) The FRITIOF Ijj] 
Monte Carlo model, based on superposition of nucleon-nucleon collisions and the Lund 
string fragmentation model, cannot describe the behavior seen in A^wtI/A™ even when 
an arbitrary rf mass modification is considered. On the other hand, hadronic cascade based 
UrQMD 16j, as well as the quark coalescence model ALCOR [17] and the thermal reso- 



nance production models of refs. 



18 



201 ] . provide a successful fit in a certain range of the 



in-medium vf masses. The main difference between the thermal models that we utilized is 
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FIG. 1. Monte Carlo simulations of A*(mT)/A™ 
The fits utilizing resonance multiplicities of refs. 
are compared to the no-mass-drop scenario of ref. 



18, 



compared with PHENIX and STAR data. 



19( (solid and dashed lines, respectively) 



JJ] (dotted line). AMPT 2.11 A(m T )/A r 



simulation with string melting 14j] report a non-thermal scenario without rj' mass modification, 
and its comparison to the data yields x 1 /ndf = 102/13, corresponding to CL = 6.8 x 10 -16 . 



in those resonance multiplicities that are not yet measured well: ref. 18 1 predicts a factor 



of 1.6 more 7/-s and a factor of 3 more rj'-s than the models of ref. 
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resonance fractions of these models are detailed in Table V of ref. [10J. Resonance decays, 



2l|. 



20| | . The relevant 



including decay chains, were simulated with JETSET 7.4 

Based on extensive Monte-Carlo simulations, x 2 of the fits to the data of Fig. [T] was 
computed as a function of m*, and B~ l for each resonance model and each fixed value 
of model parameters of a, T con d, T-po and (mt)- The best values for the in-medium mass 



of rf mesons are in, or slightly below, the range y|(2m|- + ml) < m*, < a/ 2m\ — m\ 
predicted in ref. [2J, while all are above the lower limit of m*, > y^3m n given by ref. 22]. 
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FIG. 2. Standard deviation contours on the (B , m*,) plain, obtained from A„(mT)/Af ax of 
Monte Carlo simulations based on particle multiplicities of refs. fljl . each fitted simultane- 
ously to the PHENIX and STAR combined dataset. The region between the horizontal solid lines 
indicates the range predicted in ref. [2j. The dotted horizontal line stands for Weinberg's lower 
limit |2g|. 



The A*(mT)/A™ ax simulations for the best fits of two characteristic models are compared 
to the no- mass-drop scenario on Fig. HJ while the 1, 2 and 3-cx parameter boundaries are 
indicated in Fig. [3 Those models that describe both PHENIX and STAR A*(m T )/A™ ax data 
in a statistically acceptable manner with the assumption of a sufficiently large in-medium 
Tj' mass reduction are all used for the estimation of systematics. The key parameters of the 
best fits are listed in Table HI 

Results: We have used different input models and setups to map the parameter space 
for a twofold goal, i) We excluded certain regions where a statistically acceptable fit to 
the data is not achievable, thus we can give a lower limit on the rj' mass modification. At 

6 



~ i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r~ 

no mass drop 

Letessier et al. 

— Kaneta et al. 
--- UrQMD 

total uncertainty 



-rn — i — | — rn — i — | — m — i — | — i — m — | — m — i — | — i — m — | — i — m — | — i — r-r 



from RHIC 
Au+Au 200 GeV 
central BEC data 



10' 




I 1 1 ■ I 1 ■■ I 1 ■■ I 1 ■■ I 1 ■■ I ■ ■■ I 




0.2 0.4 0.6 0.8 1 1.2 1.4 

m T -m T) [GeV] 



0.2 0.4 0.6 0.8 1 1.2 1.4 

m T -m , [GeV] 



FIG. 3. a) Reconstructed mx spectrum of the r\ and b) the rj' mesons. The dotted line indicates 
the scenario without an in-medium rf mass reduction, while the dashed, dot-dashed and solid lines 
show the enhancement required to describe the dip in the low ijit region of A* corresponding to 
the resonance multiplicities of refs. Q, Q, [l9| respectively. The shaded regions indicate the total 
uncertainty of the reconstruction. Normalization was carried out with respect to the r/ multiplicity 
of the model described in ref. jl8( |. 



Resonance m*, x 2 (CL %) 5-cr limit 

f 1 f 

model (MeV) nci/=ll * " m*, (MeV) 

< 700 

< 730 

< 570 

< 570 

< 660 

TABLE I. Most probable fits of m*, for different resonance multiplicity models with the corre- 
sponding integrated enhancement factors /„/ and of the r( and rj spectra respectively. The 
errors on the m*, values represent the 1-a boundaries of the fits. The 5-cr limits of maximum 
in-medium masses including systematics are also shown. The fitted inverse slope parameters are 
42 < B- 1 < 86 for each model. 



ALCOR [17] 490+^ 20.2 (4.29) 43.4 5.25 

Kaneta [18] 530+^ 22.8 (4.12) 25.6 3.48 

Letessier [19] 340+^ 18.9 (6.35) 67.6 4.75 

Stachel [20J 340+^ 18.8 (6.38) 67.6 4.97 

UrQMD [16] 400+^° 19.0 (6.14) 45.0 7.49 
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the 99.9 % confidence level, corresponding to a more than 5-a effect, at least 200 MeV in- 
medium decrease of the mass of the t/(958) meson was needed to describe both STAR 0-5 
% central and PHENIX 0-30% central Au+Au data on A*(m T )/Af ax in = 200 GeV 

Au+Au collisions at RHIC, in the considered model class, ii) We have determined the best 
values and errors of the fitted m*, and B~ l parameters. The best simultaneous description of 
PHENIX 5| and STAR [(J relative intercept parameter data is achieved with an rf mass that 
is dramatically reduced in the medium created in central Au+Au collisions at RHIC from 
its vacuum value of 958 MeV to 340 + f n ± 45 MeV. The first error here is the statistical 

—60 —140 

one determined by the l-cr boundaries of the fit. The second error is from the choice of 
the resonance model and the parameters (a, T con d, Tpo arid (i*t)) of the simulation. The 
third error is the systematics resulting from slightly different PHENIX and STAR centrality 
ranges, particle identification and acceptance cuts. These effects have been estimated with 
Monte-Carlo simulations, detailed in ref. jloT ]. not to exceed 9.8%, 7% and 3% respectively. 
The main source of systematic errors is the choice of the resonance models. This is due to the 
unknown initial 7/ multiplicity, hence models like ref. [l8] with larger initial 7/ abundances 
require smaller in- medium 7/ mass modification, as compared to the models of ref. {19I I20I ] . 

In addition to the characterization of the in-medium 7/ mass modification, the transverse 
momentum spectra of the r\ and 7/ mesons have also been determined. These spectra may 
serve as controls and provide motivation for future measurements as well as an input for 
theoretical calculations that may go well beyond the scope of the present manuscript. The 
low transverse momentum enhancement of the 7/ and 77 spectra corresponding to the best 
fits is shown in Fig. [3j Let us note that the enhancement of the rj production affects the 
Pt < 1 GeV region only. Our results do not modify the agreement of resonance models with 
the measured 77 spectrum in the Pt > 2 GeV region 23]. 



Discussion: Detailed analysis of the STAR and PHENIX A*(77?/r)/A™ ax dataset recorded 



at 7.7, 9.2, 11.5, 39 and 62.4 GeV during 2010 has just been started [24j, marking the 
beginning of the RHIC energy scan program. 



At present, detailed data are available from the NA44 collaboration at a/^nn = 19.4 GeV 25] 
as well as from the STAR collaboration at ^/snn = 62.4 and 200 GeV Cu+Cu and Au+Au 
collisions, the latter at different centrality classes within the 0%-80% range [8j. The NA44 
data does not feature an 77' mass drop effect. A positive sign of the 7/ mass modification 
is apparent in each case of the STAR datasets, indicating that the mass modification effect 
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is nearly at maximum in a /snn = 200 GeV Au+Au collisions and reduces with decreasing 
centrality, colliding energy and system size. We have estimated the magnitude of the system 
size and energy dependence between 62.4 GeV Cu+Cu and 200 GeV Au+Au collisions to be 
not larger than 15%, which is substantially less than the dominant systematic error coming 
from the choice of the resonance model. 

Yet to be discovered high mass resonances might also lead to an enhancement of the soft 
pion production. Indeed, such an alternative scenario can successfully explain the energy 



dependence of the K + /ti + and K~/7r~ ratios in relativistic heavy ion collisions [26J]. We 
have numerically tested the stability of our results for the enhancement of the pion halo 
coming from large mass long-lived resonances by switching on and off the contribution of 
the 0(1020) meson to estimate an upper limit of 2% on the possible effect of exotic high 
mass resonances. 

The dilepton spectrum has been measured recently in minimum bias Au+Au collisions at 
v /snn = 200 GeV, and a large enhancement was observed in the low invariant mass region 

mass enhancement of the and r] production results in 
dilepton enhancement just in this kinematic range [2|. Estimations using the enhancement 
factors in Table U indicate that the observed in-medium rf mass drop is indeed a promising 
candidate to explain this dilepton excess. 

PHENIX recently reported a two-component transverse momentum spectrum in dilepton 



channel direct photon measurements 27j], which provides an additional testing possibility to 
constrain the two component structure of the rj' and rj spectra reconstructed here. 

In summary, we report on a statistically significant, indirect observation of an in-medium 
mass modification of the rf mesons in ^/snn = 200 GeV Au+Au collisions at RHIC. A similar 
search for in-medium rf mass modification provided negative result in S+Pb reactions at 

n 

CERN SPS energies [7J. Further, detailed studies of the excitation function, the centrality 
and system size dependence of the A*(m T )/A™ ax could provide important additional details 
about the onset and saturation of the partial ?7a(1) symmetry restoration in hot and dense 
hadronic matter. Studies of the low-mass dilepton spectrum and measurements of other 
decay channels of the rj' meson may shed more light on the reported magnitude of the low 
Pt rf enhancement and the related Ua(1) symmetry restoration in high energy heavy ion 
collisions. 
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